The NaCl (100) Surface: Why Does it Not Melt? * 
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The high temperature surface properties of alkali halide crystals are very unusual. Through molec- 
ular dynamics simulations based on Tosi-Fumi potentials, we predict that crystalline NaCl (100) 
should remain stable without any precursor signals of melting up to and even above the bulk melt- 
ing point T m . In a metastable state, it should even be possible to overheat NaCl (100) by at least 
50 K. The reasons leading to this lack of surface self-wetting are investigated. We will briefly dis- 
cuss the results of calculations of the solid-vapor and liquid-vapor interface free energies, showing 
that the former is unusually low and the latter unusually high, and explaining why. Due to that 
the mutual interaction among solid-liquid and liquid-vapor interfaces, otherwise unknown, must be 
strongly attractive at short distance, leading to the collapse of any liquid film attempting to nucleate 
at the solid surface. This scenario naturally explains the large incomplete wetting angle of a drop 
of melt on NaCl (100). 
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I. INTRODUCTION 

Interest is increasing toward adhesion, the structure and physics of solid- 
liquid interfaces, and the structure of liquid surfaces, particularly of com- 
plex and molecular systems. In order to gain more insight into these prob- 
lems, there is a strong need for good case studies, to use as well-understood 
starting points. 

One of the easiest starting points may be to study the relationship and 
contact of a liquid with its own solid, a clear situation where there will be 
no ambiguity of physical description, no uncertainty in chemical composi- 
tion, no segregation phenomena, all of them complications present in the 
study of adhesion between different substances. 

Adhesion of a liquid onto the surface of its own solid usually materializes 
spontaneously with temperature. Most solid surfaces are known to wet 
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themselves spontaneously with an atomically thin film of melt, when their 
temperature T is brought close enough to the melting point T m of the 
bulk solid. The phenomenon whereby the thickness l(T) of the liquid film 
diverges continuously and critically as T — > T m is commonly referred to 
as (complete) surface melting0,Q- 

There are actually a number of exceptions to this behavior. Some solid 
surfaces in particular remain fully crystalline as T — ► T m . This sur- 
face non-melting phenomenon, originally discovered in molecular dynam- 
ics simulations of Au(lll)[!| and independently observed experimentally 
in Pb(lll)]4j, is known for the close-packed face of other metals too, such 
asAl(lll)0|. 

Here we are concerned with the surface of alkali halides, crystals well 
known for their unusually stable neutral (100) faces. Addressing a long 
time ago the NaCl (100) surface, bubble experiments by Mutaftschiev and 
coworkers revealed incomplete wetting of the solid surfaces by their own 
melt moreover with an extraordinarily large partial wetting angle of 
48°. This kind of incomplete wetting, as is physically clear, and as was 
demonstrated on metals surfaces 0, 0, is associated with non- melting of 
the crystal surface. 

In this paper we will review our recent theory work, where we showed 
by direct simulation the surface non-melting of NaCl (100) |3, Il0| . The 
reasons leading to this kind of surface non-melting are investigated. First, 
we will show the solid surface free energy calculated by thermodynamic 
integration, and see that at high temperature it drops due to a larger 
anharmonicity than in most other solid surfaces. Next, we will examine 
the surface tension of the liquid NaCl surface, we will find unusually large, 
owing to a surface entropy deficit, connected with the local surface short 
range molecular order. The solid-liquid interface free energy will finally be 
argued to be large, due to a 26% density difference. We will also discuss 
qualitatively - were hypothetically the solid-vapor interface to split into a 
pair of solid-liquid and liquid-vapor interfaces, with a thickness I of liquid 
in between - the interaction free energy V(l) expected between them. This 
interaction is here strongly attractive at short range, leading to the collapse 
of any liquid film attempting to nucleate at the solid-vapor interface, and 
causing surface non-melting. 

All quantitative results reviewed here were derived by means of molec- 
ular dynamics simulations, carried out extensively for NaCl (100) slabs. 
These simulations, as will be detailed below, are entirely based on inter- 
atomic potentials that Mario Tosi refined and published, together with 
Fausto Fumi, just over 40 years ago[12|. It is a fitting tribute to Mario's 
scientific perseverance, thoroughness, and general dependability, that these 
potentials still turn out to be so incredibly accurate, even well outside the 
range of temperatures for which they were constructed and tested, such a 
long time ago. 
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II. SIMULATIONS WITH TOSI-FUMI POTENTIALS 

The high temperature properties of solid NaCl bulk and NaCl (100) slabs 
were studied by classical molecular dynamics (MD) simulations. NaCl 
was described with the potential of Tosi and Fumi[12| who accurately 
parametrized a Born-Mayer-Huggins form. The Coulomb long-range in- 
teractions were treated by the standard three dimensional (3D) Ewald 
method applied to a geometry consisting of infinitely repeated identical 
crystal slabs. Simulated systems generally comprised about 2000-5000 
molecular units, with a time-step of 1 fs, and a typical simulation time 
of 200 ps. Long-range forces severely limit size and time in these simula- 
tions by comparison with the order-of-magnitude larger sizes and longer 
times typically affordable for systems with short range forces [3, LL3J- We 
took explicit care to ensure that all our results are not crucially affected 
by small sizes, and that full equilibration was achieved in all cases. We 
checked that 80 A of vacuum between repeated slabs are sufficient to pre- 
vent the interaction of a liquid slab with its own replicas 10]. Calculations 
were done at constant cell size and with periodic boundary conditions. 
Thermal expansion was taken care of by readjusting the (x, y) size of the 
cell at each temperature so as to cancel the (x, y) stress in the bulk solid. 
The theoretical thermal expansion was 4.05x 10~ 5 K _1 , compared with 
3.83xlO _5 K~ 1 in experiment. The theoretical bulk melting temperature 
Tm of NaCl was calculated by two phase coexistence to be 1066 ± 20 K. 
Remarkably, this Tosi-Fumi melting temperature is extremely close to the 
experimental melting temperature of 1073.8 K. The volume expansion at 
melting is about (27 ± 2)%, also in excellent agreement with the experi- 
mental value of 26%. 

Tab. d lists some of the calculated thermodynamical quantities at high 
temperature, close to the melting point. These results provide an indepen- 
dent confirmation of the outstanding quality of the Tosi-Fumi description 
of thermodynamics of NaCl, even at very high temperatures and as we 
shall see even at surfaces, where it was by no way guaranteed. 

III. SURFACE FREE ENERGIES AND NON-MELTING OF 

NACL(IOO) 

High temperature simulations of crystalline NaCl(100) slabs directly 
showed the full stability of the dry, solid surface up to Tm- Moreover, 
a well pronounced metastability of the slab solid faces above the melting 
point Tm indicated a clear surface non-melting behavior. We found that a 
much higher ("surface spinodal") temperature Tg ~ Tm + 150 K needs to 
be reached before the crystalline NaCl (100) surface spontaneously melts. 
Ts — Tm thus represents the maximum ideal overheating that a defect 
free NaCl (100) surface can theoretically sustain without becoming spon- 
taneously unstable against melting. This large overheating is quite similar 
m mag nitude to that predicted for, e.g., Au(lll) and Al (Hl)H[ll. 
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Simulation Experiment 



T M (K) 


1066±20 


1074 


AV 


27% 


26% 


L (eV/molecule) 


0.29 


0.29 


AS m (An) 


6.32 


6.38 


dP/dT (kbar/K) 


0.0311 


0.0357 


RMDS (ave.) (A) 


0.60 


0.49 


5 


20-24% 


17-20% [llj 



TABLE I: High temperature properties of NaCl. TmIs the melting temperature; 
AV is the volume jump at the melting point; L is the latent heat of melting; AS m 
is the entropy variation at the melting point; dP/dT is the resulting Clausius- 
Clapeyron ratio at the melting point. RMSD is the averaged root mean square 
displacement of atoms at the melting point; S is the RMSD over the Na-Cl 
distance, for the Lindemann melting criterion. 



At any temperature between Tm and Tg , bulk melting can only originate 
through nucleation. Even though nucleation of the melt is in reality likely 
to proceed from a localized surface droplet or defect, it is instructive to 
consider a very idealized nucleus consisting of a uniform liquid film of 
thickness I. As a function of temperature, there will be a critical nucleation 
thickness l cl -it decreasing from oo to zero between Tm and Ts 0, El ■ The 
free energy difference per unit area between a surface with a liquid film of 
thickness I and the same surface in its full crystalline state is 

G(l) = -pXl (|- - l) + ( 7SL + 7lv - 7sv) + V(l) (1) 

The first term is the gain due to the melting of the solid at T > Tm- Here 
A is the latent heat per unit mass and p is the liquid mass density. The 
second term A 700 = (7sl + 7lv^ 7sv)> is the free energy imbalance caused 
by replacing the SV interface with the SL+LV pair of interfaces, supposed 
to be non-interacting. The last term V{1) is an interface interaction, rep- 
resenting the correction to 7sl + 7lv when the two interfaces are at close 
distance. This definition implies y(+oo) = and V^(0) = — A^^. At very 
large distance the interaction disappears. The solid-vapor crystal surface 
is instead recovered when the SL and LV interfaces collapse, and the liquid 
film disappears at I = 0. The non-melting condition 7sl + 7lv > 7sv, or 
A'yao > implies that here the interaction V(l) is mainly attractive. 

This formulation indicates three possible origins for non-melting: an 
exceptionally low free energy 7sv of the solid surface; an unusually large 
free energy 7sl of the solid-liquid interface^a relatively high surface tension 
7lv °f liquid NaCl. As detailed elsewhere |9t Ho| . all three mechanisms are 
actually relevant to NaCl(100). 

The solid- vapor interface free energy at the melting point was calculated 
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FIG. 1: The solid surface free energy of NaCl(lOO) calculated from thermody- 
namic integration (circles: effective harmonic approximation). 



through standard thermodynamic integration, using 



where E is the surface internal energy, extracted from simulation of the 
crystalline NaCl(lOO) slab and the corresponding bulk at increasing T 
from 50 K to 1200 K. The surface free energy in Fig.^shows a large drop at 
high temperatures, with an increasing deviation from an effective harmonic 
behavior above 600 K, indicating very strong surface anharmonicity in this 
regime. The main source of this anharmonicity is connected with large 
root mean square thermal fluctuations of the surface CI and Na ions above 
20% of the Na-Cl distance, largely exceeding the canonical Lindcmann 
values. 

The liquid-vapor free energy, equal to the liquid surface tension, was 
evaluated from simulations of liquid NaCl slabs via the standard Kirkwood- 
Buff formula^]. The first thing we note in the result, shown in Figgis 
that right at T m the solid and liquid surface free energies are essentially 
identical, 103±4 and 104±8mJ/m 2 respectively. This is very unusual, and 
implies directly surface non-melting, because clearly A7oo = (7sl + 7lv — 
7sv) > 0. In fact, even though we did not calculate 7sl, this interface 
free energy has no reason to be very small, owing to the large solid-liquid 
density difference. We independently estimated a lower bound for 7sl to 
be 36±6mJ/m 2 10]. 

The question that remains to be explained is therefore the physical rea- 
son why the liquid surface tension is so relatively high. The liquid surface 
density profile, in particular, is very smooth, with none of the layering 
phenomena displayed by the metal surfaces (see Fig.|3|) 




G 




FIG. 3: Surface density profile of liquid NaCl, compared with that simulated 
for liquid Ar, and for two liquid metals. Note that NaCl does not show layering 
as the metals do, and has an even smoother profile than Ar 
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An important clue is provided by surface entropy (per unit area) 
SWrf = —dj/dT. The temperature dependence of two surface free ener- 
gies of Fig.[3] shows a factor 2.6 lower surface entropy Slv of the liquid 
surface compared with that of the solid surface. This liquid surface en- 
tropy deficit (SED) strongly suggests some underlying surface short range 
order. Short range order can in turn also explain why the surface tension 
is here as high as the solid surface free energy. The surface profile indicates 
that the order is clearly not layering: so what is it instead? 

The answer we find is that charge order, already very important in 
bulk, plays a newer and enhanced role at the molecular liquid surface. If 
surface thermal fluctuations are indeed very large, we find them revealingly 
correlated. For a Na + ion that instantaneously moves e.g., out of the 
surface, there is at least one accompanying CI - , also moving out; and vice 
versa. So on one hand the large fluctuations smear the average liquid vapor 
density profile, bridging very gently between the liquid and essentially 
zero in the vapor, (Fig.|3J). On the other hand the two body correlations, 
described e.g., by the the Na-Cl pair correlation function, or by its integral, 
the ion coordination number N, do not vanish identically in the vapor, but 
tend to a typical value corresponding to the NaCl molecule, (plus in fact a 
large concentration of dimers, Na2Cl2). The lack of freedom implied by the 
incipient molecular bond explains the entropy deficit, and the consequent 
large surface tension of the molten salt surface. To confirm if this is true, 
we recalculated 7sl modifying the forces in the Kirkwood-Buff formula by 
removing all contributions from surface ions whose coordination is between 
zero and 1.3, the average NaCl vapor value, which amounts to suppress 
the molecular order at the liquid- vapor interface |9(. This construction, 
meant to provide a qualitative estimate of where would the surface tension 
drop if surface molecular order were absent, gives a surface tension of 
about 50 mJ/m 2 19( . With a surface tension this low, surface non-melting 
would in fact disappear, and wetting of the solid surface by the molten 
salt would be complete. Hence the high surface tension of liquid NaCl can 
indeed be ascribed to surface molecular short range order, ultimately due 
to charge neutrality. This result confirms an early surmise by Goodisman 
and Pastor p^. 



IV. INTERFACE INTERACTION 



In this short speculative section we further rationalize the results above 
within the phenomenological framework of Eq.(l), where besides the bare 
interface free energies just calculated, an interaction V(l) appears. We 
will not present a calculation of V(l), but simply discuss it on physical 
grounds, in the light of our new microscopic understanding gained through 
simulations, and calculations of interface free energies just reviewedQ- 
The definition of interface interaction V(l) given earlier implied ^(+00) = 
and V^(0) = — Ajoo- At very large distance the interaction disappears. 
At the opposite limit, when the SL and LV interfaces approach each other 
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and merge at very close quarters, they will eventually yield the SV interface 
upon their collapse, when the liquid film disappears altogether at I = 0. 
The non-melting condition 7sl + 7lv > 7sv, or A 700 > implies here that 
the interaction V(l) is attractive at short range. 

In non-melting metal surfaces, a source of finite-range attraction was 
described as the result of a constructive interference between two equal- 
period damped density oscillations, one entering the liquid film from the 
solid side, the other, due to surface layering, from the vacuum side. Here, 
one of the two oscillations, namely that on the vacuum side, is missing, 
because there is no layering at the molten salt surface. At large dis- 
tance the main interaction between the SL and the LV NaCl interfaces 
will essentially be due to electrostatic forces and to dispersion forces. 
The latter in particular give rise to an additional long-range interface 
interaction Vdi s (0 = Hl~ 2 which is dominant at large distances. Here 
H is the Hamaker constant, that can be estimated through the formula 
H = {ir/l2)C 6 {n s - n/)(nj - n v ) = 0.00119 eV, where C 6 = 72.5 eV A 6 
is the coefficient of the Lennard- Jones interaction between chlorine ions, 
n s , m and n v are the number of CI - ions per unit volume respectively 
in the solid, liquid and vapor phases [13 . Since the liquid density is only 
about 79% that of the solid, this constant is positive which implies a long 
range repulsion of the SL and LV interfaces. A simple estimate indicates 
however that Vdi s < H/a 2 = 0.51 mJ/m 2 , a value that makes it irrelevant 
in practice. 

Therefore we expect the effective interaction V(l) to be very weak in 
NaCl, everywhere except very close to zero range, I rs a. Here it will 
suddenly turn strongly attractive, V(0) ~ — A7oo. The physics of this 
short range attraction has already been described, because it amounts to 
the free energy gained by replacing the two costly LV and SL interfaces, 
with the single and less costly SV interface. 

V. CONCLUSIONS 

Summarizing, the NaCl (100) surface is predicted to show non-melting 
and to sustain overheating up to a theoretical maximum of about 150 K 
above the bulk melting point. The thermodynamics of surface non-melting 
in alkali halides is shown to differ from that of metal surfaces, e.g. Al (111), 
Pb(lll) or Au(lll). Unlike metals, non-melting in alkali halides is not 
connected with liquid layering, but to molecular short range order raising 
the liquid surface tension, as well as to strong anharmonicty that low- 
ers the free energy of the solid surface. It is argued moreover that the 
thermodynamical SL-LV interface interaction should consist mainly of a 
strong short-range attraction. Fresh microscopic experimental work, ab- 
sent so far, is called for to check these predictions on the high temperature 
behavior of NaCl(100) and other alkali halide surfaces. 
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